Thermodynamically stable noncomposite vortices in mesoscopic 
two-gap superconductors 



o 
o 

<N 

Ph. 
< 



c 
o 
o 

Ph! 



a 
o 
o 



> 
in 

m 
m 

o 



L.F. Chibotaru, V.H. Dao and A. Ceulemans 

Division of Quantum and Physical Chemistry, and 

INPAC - Institute for Nanoscale Physics and Chemistry, 

University of Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium 



PACS 74 . 78 . Na - Mesoscopic and nanoscale systems 

PACS 74.20.De - Phenomenological theories (two-fluid, Ginzburg-Landau, etc.) 
PACS 74 . 25 . Ha - Magnetic properties 

Abstract. - In mesoscopic two-gap superconductors with sizes of the order of the coherence length 
noncomposite vortices are found to be thermodynamically stable in a large domain of the T — H 
phase diagram. In these phases the vortex cores of one condensate are spatially separated from 
the other condensate ones, and their respective distributions can adopt distinct symmetries. The 
appearance of these vortex phases is caused by a non-negligible effect of the boundary of the 
sample on the superconducting order parameter and represents therefore a genuine mesoscopic 
effect. For low values of interband Josephson coupling vortex patterns with L\ ^ L2 can arise in 
addition to the phases with L\ — L2, where L\ and L2 are total vorticities in the two condensates. 
The calculations show that noncomposite vortices could be observed in thin mesoscopic samples 
of MgB 2 . 



Two-gap superconductivity started to attract much at- 
tention in connection with the discovery of the MgE>2 su- 
perconductor, for which a clear-cut evidence for the ex- 
istence of two gaps was obtained [1]. Superconductors 
of this kind show new qualitative effects with respect to 
conventional ones. First, the presence of two or more elec- 
tronic bands at the Fermi level always enhances the super- 
conductivity as compared to the effect from any individual 
band. Second, a new phenomenon, the fractionalization 
of the magnetic flux associated to individual vortices in 
massive two-gap superconductors is predicted [2,3]. The 
condition for this fractionalization is the inequality of the 
winding numbers of the vortices in the two condensates 
(Li 7^ L2) having a common vortex core (composite vor- 
tices). Although these vortex phases have finite energies 
per unit length they never correspond to the ground state, 
i.e. are thermodynamically unstable [2,3]. On the other 
hand the energy per unit length of a vortex configura- 
tion where the vortices in each of the two bands are spa- 
tially separated (noncomposite or deconfined vortices) was 
found to be divergent [2] . These results are in line with the 
fact that only composite vortices with L\ = L2 = 1, i.e. 
usual Abrikosov vortices, have been observed in massive 
two-gap superconductors to date. 



Fractionalization of the vortex magnetic flux can also 
be observed in layered superconductors through thermal 
fluctuations [4] . By the action of the interlayer Josephson 
coupling and the magnetic field, pancake vortices (which 
individually live inside one layer) align themselves into 
vortex stacks which thread across the layers. Since each 
pancake vortex carries only « $o/Ni (where Ni is the 
number of layers) the dissociation of a vortex stack results 
in a net fraction of the flux quantum. This mechanism is 
however different from the one in a multigap superconduc- 
tor where a fractional flux vortex lives within one conden- 
sate but threads through the entire sample thickness. 

In mesoscopic superconductors, the geometry of the 
confinement of the superconducting condensate influences 
essentially their properties because the coherence length, 
£(T), and the penetration length, A(T), become of the or- 
der of the samples size (for a review, see, e.g., [5,6]). In 
particular, new vortex patterns arise as function of bound- 
ary geometry leading to cusp-like normal-superconducting 
phase boundary, T C (H) [6,7], and to the nucleation of giant 
vortices in disks [8,9] and of antivortices in regular poly- 
gons [6,10]. Another important mesoscopic effect is that 
the vortices nucleate in patterns which are quite different 
from an Abrikosov lattice and approach the latter by a se- 
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ries of phase transitions when temperature is lowered from 
T C (H) [9,11]. The symmetry of nucleated vortex patterns 
is always consistent with the geometry of the boundary of 
the sample [6], while the region of their stability against 
Abrikosov type vortices increases with diminishing size of 
the sample, persisting down to T = when the samples 
reach some critical size [11]. 

In this Letter we investigate the vortex patterns in thin 
mesoscopic two-gap superconductors. Contrary to what 
was found in the case of massive superconductors, we show 
that noncomposite vortices can arise as thermodynami- 
cally stable phases in mesoscopic samples of small enough 
size. We find the reasons for their stabilization and discuss 
the possibility of their experimental observation. 

Consider a superconducting sample of size R and thick- 
ness d in a perpendicular uniform magnetic field H. For 
small (R ~ £) and thin (d <C £ < A) samples one can 
neglect the variation of the order parameter across thick- 
ness and the distortion of the magnetic field induced by 
screening and vortex currents 1 [6,7,9]. The two compo- 
nents of the order parameter ^2 [12] are found from the 
minimization of the following 2D Ginzburg-Landau func- 
tional [13]: 
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where A is the vector potential of the applied field, 
ai = —ait and ot2 — <^20 — are the condensation 
energy parameters for the active and the passive band, 
respectively, t ~ 1 — T/Ti, and T\ is the critical tem- 
perature corresponding to the active band [13]. At the 
normal/superconducting phase transition, minimizing (1) 
(without the terms ~ |^n| 4 ) results in two linear equa- 
tions describing the nucleation of superconductivity and 
two boundary conditions: 
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where n,n'=l,2 and 2,1. The notation | n .b. means the 
expression is projected on the unit vector normal to the 
boundary. It is straightforward to show that the nucle- 
ation solution of these equations has the form ~ ^2 ~ 
(f>N 1 where 0at is the solution of the eigenvalue equation: 

2e 2 

i-ifiV A) <j>i = Xifc, 
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1 For a thin superconductor with ^>(r) — y 
variation of the vector potential due to the supercurrent (5A(i*o) « 
-4&r / (f?Vx+A)/ 2 (r)/|r - r|dS ~ (§£V X +A)/ 2 (r ) can be 

neglected if the thickness d is small enough. The same argument can 
be generalized to a two-gap superconductor. 



corresponding to Xn(H) which is the lowest of the eigen- 
values Xi(H). This eigenvalue determines the nucleation 
phase boundary T C (H) [13] and the ratio of the ampli- 
tudes of the two components of the order parameter, 
^2/^1 = 7/(<*2 + AAr/2m 2 ). Thus the distribution of the 
nucleated order parameter in the two bands of a two-gap 
superconductor is described by the same function. 

The eigenstates of (3) are then used as the basis set for 
the order parameter of the complete functional (1), 



(4) 



which yields AF as a function of the expansion coefficients 

We will consider further a disk of radius R. It is 
then convenient to introduce new "lengths" defined by 
the relations: —a\ = ft 2 /2mi£i(T) 2 , o^o = ^V^i^fo? 
-a 2 t = h 2 /2m 1 &(T) 2 , and 7 = H 2 /2m 1 ^. When 7 = 
and 0L2 > the superconductivity nucleates in the active 
band, for which the coherence length is £1. The order pa- 
rameter corresponds to one single giant vortex in the cen- 
ter of the disk [8] whose winding number increases with the 
applied field following the cusps on the nucleation phase 
boundary line (Fig. la). At lower temperatures it under- 
goes a broken-symmetry phase transition to a multivortex 
state corresponding to Abrikosov (L = 1) vortices form- 
ing a regular polygon [9]. Since 1/^ ~ T c — T, it follows 
that the temperature region where the nucleated (giant 
vortex) phase is thermodynamically stable (dark regions 
in Fig. la) scales with the size of the disk as ~ R~ 2 . 
Note that for R > 6£i new broken- symmetry phases arise, 
which correspond to off-center displacements of the cen- 
tral vortex in the case L = 1 and of the maximum of |^| 2 
distribution in the case L = [14]. 

In the case of nonzero interband Josephson coupling 
similar phase diagrams emerge (Fig. lb). We see again 
regions corresponding to giant vortex patterns (shown in 
blue in Fig. lb), with the same vortex numbers in both 
condensates, and regions with broken-symmetry vortex 
patterns (yellow and pink regions in Fig. lb) consist- 
ing of Abrikosov vortices. However a qualitatively new 
feature arises in the broken-symmetry phases, with two 
vortex patterns, corresponding to ^1 and ^2, being spa- 
tially distinct for any parameters of (1). Therefore each 
vortex in these phases is a noncomposite vortex. Moreover, 
in contrast to single band superconductors (Fig. la), we 
have now two types of multivortex phases, corresponding 
to Li = L2 (yellow regions in Fig. lb) and L\ ^ L2 (pink 
regions in Fig. lb). The corresponding order parameter's 
density plots are shown in Fig. 2. 

The temperature evolution of the vortex patterns in two 
condensates is strongly dependent on the parameter a 2 , 
while the difference in their shape increases with dimin- 
ishing 7. Indeed, as Fig. 2 shows with i?/£ 7 (0) = 0.1, 



p-2 



Thermodynamically stable noncomposite vortices in mesoscopic two-gap superconductors 



30 



20 



DC 



10 



30 



u7 20 

C\J 

DC 

10 



I I I I I I I I I I I II I I I I I I I I II I I I I 

■(a) / /© 


i i i i i i i i_ 

©1 




; ° / ri 






/ A®)/ 

/ / V J / 






q ® ® 


® 




i \\\\ 

l\r<\ II LI' '\ 






■ o -/pi 






t 1m 







3 4 | 5| ^ 6 

OM> 2-(a) 2-(c) 2-(b) 



Fig. 1: (color) Phase diagram for a thin disk of single-gap 
(a) and two-gap (b) superconductors. The phase diagram (a) 
is parameter- free [11], while the diagram (b) is for a\/a2 — 
fr/fc = mi/m 2 = 1, R/&(0) = 2^3 and i*/f 7 (0) = 0.1. 
For each phase, the vortex structure is shown schematically in 
black for composite vortices and in green (active band) and red 
(passive band) for noncomposite ones. The numbers inside the 
circles indicate the vorticities of the giant vortex phases. Verti- 
cal lines separate phases with different total vorticity L. Con- 
tinuous horizontal lines correspond to broken-symmetry phase 
transitions, which are of giant vortex - multivortex type for 
L > 1. The dashed line is the nucleation phase boundary for 
the passive band for 7 = (see the text). The vertical dotted 
lines mark the positions of the three graphics in Fig. 2. 



for small values of 7 we encounter three qualitatively dif- 
ferent situations. In the first case (a) the transition to 
the broken- symmetry phase leads to a strong separation 
of the multivortex patterns, corresponding to the active 
(green) and the passive (red) bands, the latter approach- 
ing the first one with further lowering of temperature. In 
the second case (b) the vortex pattern in the passive band 
first shrinks, resulting in its strong separation from the 
vortex pattern in the active band, and then approaches it 
again when the temperature is lowered. Finally, the case 
(c) corresponds to the transition to the broken-symmetry 
phase phase with L\ 7^ L2, followed by the transition to a 
broken-symmetry phase with L\ — L2 when the tempera- 
ture is lowered further. 

The reason for the above behaviour of the order param- 
eter can be elucidated by the following simple considera- 




Fig. 2: (color) (a) Evolution of the order parameter for three 
values of the applied magnetic flux in the phase diagram of 
Fig. lb. The background colours correspond to the giant vor- 
tex phases (blue) and the multivortex phases with L\ — L2 
(yellow) and L\ 7^ L2 (pink). The expansion coefficients ul 
and vl correspond to the ground state solutions of (3) for each 
L that contribute substantially to the expansion. The numbers 
in the corner of inserted plots gives the side of the zoom on the 
disk in units of R. The increase of the intensity of colors corre- 
sponds to the decrease of the density of the order parameters, 
so that the darkest regions denote the position of the vortex 
cores. 



tion [11]. Let's suppose for simplicity that only one eigen- 
state of (3) (~ <Pa) admixes to the nucleated order param- 
eter (~ cj)jsf) in the point of broken-symmetry transition, 
which is actually a reasonable approximation for disks [9] 
and regular polygons [11]. When 7 is small the nucle- 
ation of superconductivity takes place at — ql\ ~ A at /2m\. 
At the onset of the transition to the broken symmetry 
phase we have ^1 = un^n + ua^a (<^a is supposed to 
be of different symmetry compared to 0at, which in the 
case of disk implies different rotational quantum num- 
bers). If «2 + Aiv/2m2 > (the sufficient condition for 
this is £2(0) > £20)? the minimization of (1) gives for the 
expansion coefficients of ^2' 
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where = XiR 2 /h 2 are dimensionless eigenvalues of (3), 
only £1 is temperature dependent, and the dependence on 
R is explicit. If Vi and ui were related by the same pro- 
portionality coefficient, the order parameters ^1 and ^2 
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would have been described by essentially the same func- 
tion, in particular, all vortices arising in such a phase 
would correspond to composite ones. Such a situation 
obviously occurs in the giant vortex phase. However, as 
Eq.(5) shows, the transition to a broken-symmetry phase 
makes the coefficients of proportionality different for the 
admixed (i = A) and the nucleated (i = TV) components 
of the order parameter due to the term ~ in the de- 
nominator of (5). This term becomes important for small 
values of i?, i.e. in the mesoscopic regime, and it disap- 
pears in the macroscopic regime R — > oo. This is in line 
with the conclusion mentioned in the introduction, that 
only composite vortices can arise in thermodynamically 
stable phases of massive two-gap superconductors. 

The above analysis allows to explain the phase diagram 
in Fig. lb. For large positive ot2 + the supercon- 

ductivity in the passive band is induced by the Josephson 
coupling to the active band. Then the order parameter 
and the resulting vortex pattern in the passive band fol- 
low closely that of the active band, according to Eq. (5). 
However when this quantity is small and can turn to neg- 
ative at some temperature T* (shown by dashed line in 
Fig. lb) then for T < T* the nucleated component <j)N 
in the passive band begins to grow much faster than in 
the previous case due to the intrinsic superconductivity 
which now exists in the passive band. Then two quali- 
tative situations can occur. If at T = T* the supercon- 
ductor is still in the giant vortex phase, the subsequent 
transition to the multivortex phase in the active band will 
not induce a similar order parameter in the passive band. 
The latter will remain almost ~ c/)n and consequently a 
much smaller amplitude of the splitting of the giant vor- 
tex into Abrikosov vortices is expected (Fig. 2a). If at 
T = T* the system is already in the multivortex state, 
then the subsequent lowering of temperature will induce 
the shrinkage of the vortex pattern in the passive band 
due to a faster stabilization of the nucleated component 
in ^2- This is exactly what is seen in Fig. 2b. In both 
these situations a strong spatial separation of the vortex 
patterns is achieved. The existence of the vortex patterns 
with L\ 7^ L2 is merely due to the fact that for 7 = 
the phase diagram for the passive band is shifted upwards 
relative to the active one (the nucleation of superconduc- 
tivity takes place at lower temperatures), which creates 
overlap regions where the two condensates have different 
voticities. If 7 is smaller than the difference of free ener- 
gies corresponding to L\ and L2 in the passive band then 
vortex patterns with different vorticities in the conden- 
sates will be stabilized also for 7 7^ 0, but will disappear 
from the phase diagram when i?/£ 7 exceeds some critical 
value. The existence of the overlap regions with L\ ^ L2 
is possible due to the fact that the lines separating do- 
mains of different vorticity on the diagram on Fig. 1 are 
not vertical. As a consequence, lowering the temperature 
at fixed applied flux can result in several phase transitions 
with change of vorticity in the condensates, accompanied 
by jumps of the order parameter, i.e. of the coefficients Ui 
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Fig. 3: (color) The phase diagram for a thin disk of MgE>2 
with the radius R=70 nm and the parameters taken from Ref. 
[13, 16]. The upper panels are zooms on the distribution of 
the order parameter in broken-symmetry (multivortex) phases. 
The conventions follow Figs. 1 and 2. 

and Vi in Eq. (4) (Fig. 2c). 

Finally, we calculate the T — H phase diagram for a 
thin disk of MgE>2 with the anisotropy axis c perpendic- 
ular to the plane of the sample 2 (Fig. 3). Formation of 
noncomposite vortices is favored by a weak interband cou- 
pling 7 and a small mass ratio 7x12/ m\. From current 
available data, the smallest 7 is obtained for a sample 
of Boron isotope Mg 10 E>2 with a value reduced by 30% 
from usual Mg 11 B2 [15], while m^/rria- < 0.1 has been 
observed in irradiated samples [16] (where the introduced 
disorder should change effective masses but not the cou- 
pling constants). Modifying consequently the mass ra- 
tio in the Ginzburg-Landau parameters estimated previ- 
ously for pristine MgB2 [13], we use 0,2/0*1 = P2/P1 = 1-5, 
m 2 /m 1 = 0.07, R/£i(0) = 10, R/&o = 8.1, i?/f 7 = 6.3 
and T c = 39K for Fig. 3 [17]. We can see several regions 
on the phase diagram where the two vortex patterns are 
well separated in space, some of them being shown in Figs. 
3 a-c. The radius of the disk was taken R = 70 nm but 
similar separations of the vortex patterns (3-5 nm) were 
found for a wide range of radii: R=30 - 120 nm. We have 
observed that other sets of parameters are not critical for 
the separation of the two vortex patterns except for the 
ratio m2 /mi', a smaller ratio leads to a higher separation. 

The existence of noncomposite vortex patterns in two- 
gap superconductors can be experimentally verified by 
combining different local probe techniques [18,19]. The 
discussed effects are not related to the symmetry of the 
samples but arise due to non-negligible influence of the 
boundary of the sample on the two superconducting con- 

2 This is the way in which thin MgB2 films actually grow on the 
substrate. 



p-4 



Thermodynamically stable noncomposite vortices in mesoscopic two-gap superconductors 



densates. In this sense the emergence of noncomposite 
vortices found in the present study represents a true meso- 
scopic effect, not observable in infinite superconductors. 
We note in this connection that noncomposite vortices 
could also be stabilized in ultracold atomic Fermi gases 
in optical lattices or in single traps where the BCS to 
Bose-Einstein condensation transition takes place. These 
finite-size systems possess tunable interaction parameters 
which can give rise to multiband superconductivity [20] 
and are thus the subject of new promising experiments. 

* * * 
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scattering is argued to explain the discrepancy, the ori- 
gin of the interband- coupling variation remains an open 
question. 
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